On the validity and breakdown of the Onsager symmetry in mesoscopic conductors 

interacting with environments 
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We investigate magnetic-field asymmetries in the linear transport of a mesoscopic conductor 
interacting with its environment. Interestingly, we find that the interaction between the two systems 
causes an asymmetry only when the environment is out of equilibrium. We elucidate our general 
result with the help of a quantum dot capacitively coupled to a quantum Hall conductor and discuss 
the asymmetry dependence on the environment bias and induced dephasing. 
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More than two decades ago it became clear [l[ that the 
Onsager-Casimir symmetry relations 0, Q are crucial to 
understand the transport properties of mesoscopic con- 
ductors. These symmetries are fundamentally a conse- 
quence of microreversibility of the scattering matrix that 
describes the conductance of a phase-coherent conduc- 
tor, dictating that the two-terminal linear conductance 
G must be symmetric under reversal of the external mag- 
netic field B. An interesting consequence is that the 
phase of the Aharonov-Bohm conductance oscillations of 
a solid-state interferometer with a quantum dot embed- 
ded in one of its arms can take the values or tt only, 
thus leading to phase rigidity [1, Q . 

The scattering approach to mesoscopic transport as- 
sumes that the mesoscopic conductor preserves the elec- 
tron quantum phase while inelastic processes giving rise 
to irreversibility take place only in the reservoirs that feed 
and draw the current. Therefore, close to equilibrium G 
is a function of the transmission T evaluated at the Fermi 
energy Ep common to all terminals and the Onsager sym- 
metry implies T{B) = T{—B). However, the system in- 
evitably interacts with the environment which may give 
rise to irreversible processes. Theoretically, the Onsager 
symmetry has been proved to be valid for an isolated con- 
ductor. Therefore, it is an interesting question whether 
a conductor interacting with its environment still fulfills 
the symmetry. In this Letter, we predict that this inter- 
action leads, in fact, to magnetic-field asymmetries (or, 
briefly, magneto-asymmetries) when the environment is 
driven out of equilibrium. In addition, our theory con- 
firms why magnetic-field symmetries are preserved in 
previous experiments on two-terminal conductors even 
if they cannot avoid interactions with its environment. 

Recent works 0] have shown that magneto- 
asymmetries arise in nonlinear mesoscopic transport, a 
fact which has been observed experimentally [1, 0, [3, 
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111 . Il2l Il3l |. In contrast, here we address the magneto- 
asymmetry of the linear mesoscopic conductance when 
the environment is out of equilibrium. 

Consider the model system sketched in Fig. 1. System 
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FIG. 1: Schematic representation of the system under consid- 
eration. System C is a conductor capacitively coupled to an 
environment (system D). 



C is a mesoscopic conductor coupled to reservoirs L and 
R. The environment, denoted as system D, is modeled 
as a second conductor in close proximity with system C. 
There exists a Coulomb interaction coupling conductor 
and environment electrons but no particle exchange is 
permitted between the two subsystems. Experimentally, 
the environment can be a quantum point contact (QPC), 
a quantum Hall bar or any other system whose electron 
states depend on the electronic trajectory across the con- 
ductor. The environment can be driven out of equilib- 
rium with applying a bias between reservoirs X and Y . 
As a consequence, we must consider scattering of two 
particles described by the following (asymptotic) states: 
|1) = \L)c®\X)d. |2) - \L)c®\Y)d. |3) - \R)c®\X)d 
and |4) = \R)c ® \Y)d- \ - ■ ■)c and | ■ ■ • )d represent the 
electron states at system G and D, respectively. 

Suppose that an electron from lead L is injected into 
the conductor. Before scattering, the initial density ma- 
trix Pin is given by 



Pin - {\L){L\)c ® {nx\X){X\ + nY\Y){Y\)i 



(1) 



where ux and ny are the transport electrons in sys- 
tem D that originate from leads X and Y , respectively. 
Normalization condition, Trpin — 1, gives the constraint 
nx + ny = 1. At equilibrium, one has nx = ny = 1/2 



whereas out of equilibrium we rewrite pin as 



Pn 



6p)\l){l\ + \{l-5p)\2){2\, (2) 



where Sp = nx — ny represents a nonequiHbrium pa- 
rameter. At small bias, Sp is proportional to the bias 
voltage while very far from equilibrium Sp = I. The 
limit of Sp = 1 is discussed in Ref . fl^ . Upon scatter- 
ing, the output density matrix is pout = SpinS'', 
where S is the two-particle scattering matrix with ele- 



ments 



(ij = l,-.-,4). Note that S de- 



scribes scattering between two electrons in different con- 
ductors [31. We emphasize that S cannot, in general, be 
written as the product of two single-particle scattering 
matrices and that this nonseparability is prec isely due to 
the interaction between the two systems 17 1. 

The transition operator, T, gives the transmission 
probability for an electron in system C, T = Tr (pout2^)- 
It reads r= i\R){R\)c (S Id = |3)(3| -I- 14}(4| , where /r, 
is the unit operator in system D. Then we find 

T = 1 - nx{\Snf + \S21f) - nyilSuf + \S22f) , (3) 

where unitarity of 5* has been used. 

We introduce the magneto-asymmetry factor a = 
T{B)-T[-B). Microreversibility implies that Sij{B) = 
Sji{~B). Hence, we infer from Eq. ^ that 



= Sp{\S 



121 



(4) 



Quite generally, one has 15*121^ ^ l'5'2i|^ and then the 
transmission through the conductor is clearly not sym- 
metric under reversal of the magnetic field. Notably, the 
asymmetry is proportional to Sp, thus Eq. (|3|) predicting 
that at small bias the magneto-asymmetry in the con- 
ductor grows linearly with the voltage applied in system 
D. Below we perform numerical simulations in a realistic 
system that confirms this prediction. 

Our argument can be easily extended to the case where 
the nonequilibrium situation of system D includes more 
than two leads. We note that Eq. ^ is valid when 
electrons are interacting within each subsystem and also 
when the electron at system C interacts with more than 
one electron at D but the problem then becomes involved 
because one should resort to a multiple-particle scatter- 
ing matrix. Interestingly, the magneto-asymmetry van- 
ishes when the environment is in equilibrium (Sp = 0), 
and it does not depend on a specific model for the envi- 
ronment. This explains why all linear-transport exper- 
iments satisfy the Onsager symmetry even though the 
conductors cannot be isolated from uncontrolled interac- 
tions with their environments. 

It is clear from Eq. (|4]) that the magneto-asymmetry 
is nonzero only to the extent that S is nonseparable, i.e., 
|'S'i2p = 1521 P if S is given by the product S = S'c So 
where Sc and Sd are the single-particle scattering ma- 
trices of system C and D, respectively. Unitarity of Sd 




FIG. 2: Sketch of a quantum dot capacitively coupled with 
the top edge of a quantum Hall bar for B > 0. The edge state 
is transmitted at the QPC with probability T^. For B < 
the arrow directions are reversed. 



is necessary in deriving this relation. In Ref. [18[, S is 

expressed in terms of the scattering matrix of the un- 
coupled systems to leading order in the interaction cou- 
pling strength and the correction term of the transmis- 
sion probability of the first system is found to depend on 
the injectivity of lead X of the second system. But the 
injectivity alone is not invariant under field reversal [^. 
Therefore, T need not be an even function of B due to 
interaction between the two systems. 

We focus on the zero temperature case for simplicity 
and assume that electrochemical potentials of leads X 
and Y are Ep + cVd/^ and Ep — eVo/S with Vd the 
bias of system Z?. In a two-dimensional conductor, nx oc 
Ep + cVd/^ and ny oc Ep — cVd/^. Using the constraint 
nx + ny = 1, we obtain the nonequilibrium parameter 
Sp = nx — ny — eVD/2Ep. 

Let us now illustrate the general result given by Eq. ^ 
with an instructive example as depicted in Fig. [2l 
We consider resonant tunneling through a quantum dot 
which is capacitively coupled to the top edge of a quan- 
tum Hall conductor, which works as the controllable en- 
vironment. We assume that the filling factor is = 1 and 
thus current is carried by two edge states along opposite 
sides of the sample. 

In addition, we consider a QPC constriction that par- 
titions the current injected from lead X with probability 
Td- For the moment, consider the case Td = I. Qual- 
itatively, it is clear that the current traversing the dot, 
even close to equilibrium, is not an even function of B. 
This follows from the fact that the potential of the edge 
state Ue is in equilibrium with the electrochemical po- 
tential of the injecting lead. Then, Ue = Vx for i? > 
and Ue = Vy for _B < 0. Since the conductance through 
the dot depends on the local potential at the dot Ud, and 
Ud, in turn, depends on Ue via the capacitive coupling, 
we must have, quite generally, G{B) ^ G{~B) flit . 

We investigate the case in which Ud reacts continu- 
ously to a change in Ue and treat interactions in a mean- 
field way. We describe the scattering through the dot 
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FIG. 3: (Color online) Magnetic-field asymmetry of the linear 
conductance through the dot shown in Fig. [2] for Ep = eo = 
0, r_L = Fij = r/2 — 0.1 and ef/oq = F as a function of 
the voltage bias in the Hall bar Vh- Upper inset: Linear 
conductance at different polarizations of the magnetic field. 
Lower inset: Self-consistently calculated screening potential. 
In the insets we take Cn ~ 0.1 and eVn = F. 



with a Breit-Wigner resonance with level position Eq and 
broadening F = -t-Fj^. F^ (F^j) denotes the resonance 
broadening contribution from the left (right) lead. To 
be definite, we take Eq and F invariant under B reversal, 
which is true at equilibrium. To separate equilibrium 
and nonequilibrium contributions, we consider the po- 
tential Ud = Ucq + AU. Then, away from equilibrium 
the screening potential AU follows from the charge neu- 
trality condition which establishes that the net charge 
— 'Zncq — Icq must bc equal to the polarization charge 
permitted by electrostatics: 



SqiVL,VR, Ud) = C{U,q + AU- Ue) 



(5) 



The charge injected from the left and right elec- 
trodes, g„eq = /f'^^'''^ Dl{E - eC/eq - eAU)dE + 
J^^'^'^^"' Dji{E-eUcq-eAU)dE, is given by the injectiv- 
ities Da{E) = (eV27r)F„/|A(£:)|2 where a = [L,R) 
and A{E) = E ~ Eq + iT/2. The equilibrium charge, 
Qoq = Dd{E — eUcq)dE, depends on the total den- 
sity of states, Dd = Dl + Dji. In Eq. ([5]), C is the 
geometrical capacitance between the edge state and the 
dot. If the density of states of the edge state is much 
larger than C/e^, we simply have Ue = Vx for S > 0. In 
the general case, we find 



CUd + e^DeVx 
C + e^De 



(6) 



for i? > 0. In the last equation, Vx should be replaced 
with W for S < 0. Upon inserting Eq. ^ in Eq. jS]), 
we obtain 5q = C^(J7eq + AU - Vx) where C^^ = C^^ + 
(e^Z^e)"^ is the electrochemical capacitance. 



Equation ^ is to be solved self-consistently. Once we 
numerically find AC/, we can assess the current. 



h JE^+eV^ \^{E - eC/eq - 



eAU)\- 



:dE, (7) 



and the linear conductance G = dl /dV\v=o with V = 
Vl — Vr. Upper inset of Fig. [3] shows results for G 
as a function of the equilibrium level position eUcq for 
a voltage bias ~Vx = Vy = Vh/2 = F/2e apphed 
in the Hall bar. We observe that G differs for oppo- 
site B orientations. The reason for the asymmetry is 
uniquely due to the asymmetry of the potential AU (see 
lower inset of Fig. [3]) arising from the asymmetry of the 
Hall bar injectivity. Thus, we present in Fig. [3] calcu- 
lations of the dimensionless magneto-asymmetry factor 
a = [GiB) - G{-B)]/[G[B) + G{-B)] as a function of 
Vh for eUcq = F and various G^. These results are cen- 
tral to our discussion. The magneto-asymmetry is larger 
for larger G^ since the interaction coupling to the edge 
state is stronger. Of course, in the limit C^^ — > (which 
amounts to C ^ 0), the magneto-asymmetry vanishes, 
fulfilling the Onsager symmetry relation. It also vanishes 
in the limit of an equilibrium environment {Vh 0). 
Moreover, a is linear with Vh for small Vh, in excellent 
agreement with Eq. ([4]) 21]. Saturation in a takes place 
for large Vh and C^^, for which the precise form of the 
local density of states starts to play a role. 

Let us now consider the case where the dot is coupled 
to a partitioned edge state. There is a probability Td 
{Rd = 1 — Td) that the edge state is transmitted (re- 
flected) from lead X to Y {X) through the QPC. For 
nonzero G, we find 



UeiB > 0) 



e^D,iTdVx + RdVY) + GUd 



G 



(8) 



For B < Q one replaces Vx with Vy in Eq. As a 
result, Ue is B-asymmetric for nonzero Td (for Td — Q 
the symmetry is restored) and depends in a generic way 
on Td only for i? > 0. 

More interestingly, this geometry gives rise to magneto- 
asymmetry of dephasing. Dephasing in the dot is caused 
by current partition in the Hall bar [H, [H, [H, HI, IH, 



271 . 128l ] . It is intimately related to the possibility of ex- 
tracting charge state information of the dot from the rela- 
tive phase shift between the transmitted and the reflected 
beam at the QPC. For B > 0, Coulomb interaction in- 
duces dephasing because only the transmitted electron 
undergoes a phase shift. For B < 0, however, dephasing 
is not induced because the dot interacts with the Hall bar 
before the electrons arrive at the QPC. Therefore, it is 
clear that the magneto-asymmetry of linear transport is 
caused not only by the asymmetry of the local potential 
of Eq. IH]) but also by the asymmetry of dephasing. 

Dephasing induces an additional broadening for B > 0: 
F ^ F -|- F^ with F^ = rjTdRdVn , where rj is a constant 



4 
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FIG. 4: (Color online) Magnetic field asymmetry as a function 
of the QPC transmission for various biases applied to the Hall 
bar Vh- Parameters are Ef = eo = 0, Fi = F_b = F/2 = 0.1, 
= 0.5, and e(/cq = F. Inset: Asymmetry versus Vh for 
different QPC transmissions. 



which depends on the details of the interaction between 
the edge state and the dot [H HI, HI, [ll, [111 (r? = 
for B < 0). Since our goal is to offer a simple picture of 



the effect, we adopt the phenomenological voltage probe 
model [29| . It assumes a fictitious voltage probe attached 
to the dot with coupling T^. The condition = deter- 
mines the potential at the probe, V^. Every carrier which 
enters the probe, is reemitted into the dot with a com- 
pletely unrelated phase, thereby giving rise to dephasing. 
Thus, we add a term /i^^"^''^* D^E - eU^q - eAU)dE 
in the left-hand side of Eq. ((5]) and a current contribu- 
tion {2e/h) Je^+^v^ TlT4,/\A{E - eC/cq - eMJ)\^dE to 
Eq. ([7]). Figure m shows the effect of current partitioning 
in a. The amount of dephasing is tuned with Td- (To 
emphasize the role of Td in the asymmetry we use a small 

= 10-3). When Td = 0, J7e = Vy independently of the 
B direction. As a result, 5 = 0. When Td increases, 
a enhances monotonically. For higher Vh, <5 becomes 
larger in agreement with Fig. [3l Finally, in the inset we 
plot a as a function of Vh for decreasing values of T^, 
which also demonstrates that a vanishes for Vh —* 0. 

In conclusion, the statement that the two-terminal lin- 
ear conductance must be symmetric under reversal of the 
magnetic field is widely accepted and has been exhaus- 
tively confirmed. However, conductors inevitably inter- 
act with the external environment. We have shown that 
a magnetic-field asymmetry appears even in the linear 
response when the environment is out of equilibrium. 
This situation can be realized with another conductor 
in close proximity applying an electric bias across it. Im- 
portantly, we predict that the asymmetry depends on 
the two-particle scattering matrix. We have examined a 
quantum dot coupled to a quantum Hall bar and found 
that the asymmetry grows with the Hall bar bias and the 
interaction coupling strength. It also leads to an asym- 



metry of dephasing when the Hall current is partitioned. 
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